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Abstract 

m 
vo 
o 

■ One-dimensional (ID) model of non-relativistic particles with inverse-square 

o , 

interaction potential known as Calogero-Sutherland Model (CSM) is shown to 
possess fractional statistics. Using the theory of Jack symmetric polynomial 

a 

i . the exact dynamical density-density correlation function and the one-particle 

Green's function (hole propagator) at any rational interaction coupling con- 
stant A = p/q are obtained and used to show clear evidences of the frac- 
^ . tional statistics. Motifs representing the eigenstates of the model are also 

constructed and used to reveal the fractional exclusion statistics (in the sense 
of Haldane's "Generalized Pauli Exclusion Principle"). This model is also 
endowed with a natural exchange statistics (ID analog of 2D braiding statis- 
tics) compatible with the exclusion statistics. 
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The quantum particles obeying fractional statistics known as anyons have been a subject 
of intense study since the discovery of the fractional quantum Hall effect and the High T c 
super-conductors [|IJ. The subject, however, is far from complete. In fact, even the simplest 
model of anyons, namely the ideal gas, has not yet been fully understood. To add to the 
confusion, there are now two seemingly non-equivalent definitions of fractional statistics. 
While the popular definition of anyons is based on the quantum phase arising from the 
exchange of particles Haldane's definition is based on so called the "Generalized 
Pauli Exclusion Principle." The main difference in the two approach is that while in the 
former the statistics is assigned to the Newtonian point particles, in the latter it is obeyed 
by the elementary excitations of condensed matter system. 

There are no fully solvable two-dimensional models where the ideas of fractional statis- 
tics can be rigorously tested. In one-dimension, however, there is such a model known as 
Calogero-Sutherland Model (CSM) 0]. I show in this Letter that the two definitions of frac- 
tional statistics can coexist in CSM without any inconsistency. Following Ref. Q], I construct 
the motifs for all the excited states and explicitly demonstrate that the quasi-particles and 
quasi-holes obey the Haldane's exclusion statistics. I also solve the exact ground state dy- 
namical density-density correlation function (DDDCF) and the hole propagator part of the 
one-particle Green's function (HPOGF) and show that the contributing intermediate states 
involve only a finite number of quasi-particles (holes) consistent with the ideal "anyon" 
gas structure of this model. In calculating the correlation function and propagator a new 
mathematical technique based on the theory of Jack symmetric orthogonal polynomial || 
is used. 

The CSM Hamiltonian, which describes a system of N non-relativistic particles interact- 
ing with inverse-square exchange, is given by 

x- E d \ i £ 2 8 A(A " 1) , (i) 

j=i & x j j<i d 2 (xj — xi) 

where h 2 /2m = 1 and d(xi — yj) is the chord distance between the ith and jth. particles on 
a ring of length L, and is equal to \(L/tt) sin(7r(xj — yj)/L)\. The dimension-less interaction 



coupling constant A is a positive real number that specifies the natural statistics of this 
model. For the special values of A = 1/2, 1, and 2, the model is related to the orthogonal, 
unitary and simplectic random matrix theory ||, respectively, and DDDCF for those values 
have previously been found by the supersymmetry technique J/J. For A = 2, the full one- 
particle Green's function has also been found [§,0. 

One of the main results of this Letter is the calculation of the dynamical density- density 
correlation function at any rational interaction parameter A = p/q. It is given by 

(0|p(x,t)p(0,0)|0) =C Ref[ (TdXi) f[ ( f dy 3 ) Q 2 F(q,p, \\{x uyj })S Qx - Et \ (2) 

i=i ^ Jo ' j=i \ Jo ' 

where Q (momentum) and E (energy) are given by 

Q = 27rp(t / x j + f / yj) (3) 
\j=i i=i / 

E = {2Kp)*(£ep{z J )+£e H (y J )), (4) 

\j=l 3=1 J 

where p = N/L, ep(x) = x(x + A) and €h(v) = Ay(l — y). Xj(ep) and 2/j(e#) are nor- 
malized momentum (energy) of the quasi-particles and the quasi-holes, respectively. The 
normalization constant C is given by 

A( 3, r™(A)r»(i/A) - / r( g -(j-i)/A) \ 2 

A{m,n,A) n^r^-ACi-i^ny^r^-O'-iVAj/iAra-Cj-iVA); ' [b) 

C= \ 2 ]] [P) A(q,p,\). (6) 
Finally, the form factor F(q,p, \\{xi, y^}) is given by 

The results of Simons, et al. at A = 1/2, 1 and 2 with an appropriate change of variables 
agree with Eq. (|2|) up to the normalization constant. Recently, based on a mapping to 2D 
QCD Minahan guessed correctly the form factor for the integer values of A and their duals 
1/A flQfl . Haldane also guessed the correct form factor for the rational A = p/q based on the 
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information given by Simons, et al., Galilean invariance, and asymptotic agreement with 
U(l) conformal field theory fTT |. 



The ground state of SCM is given by 

N 

*o=ife--^) A rK o (8) 

3>l j=l 

where Zj = exp(i27cXj/L) and Jo = — A(iV — l)/2. If a general eigenfunction with energy E 
is written as \1/ = then $ is an eigenstate of the following new effective Hamiltonian 

H' 

= £MJ 2 $ + A£^(^ - = ( 9 ) 

i i<j Z i Z j 

where e = (L/2tt) 2 (E - E Q ). An amazing coincidence happens if A = (AiV — A — 1) J2j z j9 Zj 
is added to Eq. (|9|). A complete set of linearly independent solution of the resulting equation 
is known in mathematical literature as the Jack polynomial J^ x (zi, . . . , zn) ||. The index 
k = K2, ■ ■ ■ , Kjv) is partition of non- negative integers and is essentially a set of bosonic 
quantum numbers used to label all the eigenstates of SCM up to global Galilean boosts. In 
particular, the ground state is given by the partition of zero (i.e. k± = «2 = ■ ■ ■ = «iv = 0). 
The parts k^s of partition k are ordered so that K\ > k,2 > • • • > k<n- Since the J K is a 
homogeneous symmetric polynomial of degree \k\ = J2j K ji it is also an eigenfunction of A 
with eigenvalue (AiV — A — and thus of H' with e = Y,j K j + ^(N + 1 — 2j)nj. If the 
pseudo-momenta k/s are defined by 

Lkj = 2nlj + tt(A - 1) 5>gn(^ - h), (10) 

with Ij = Kj + (N + 1 — 2j)/2, then the eigenenergy and momentum take the following 
free form: E = 7^J2j kj an d P = J2j kj. Using a method developed by Yang and Yang 
fi~2|| , Sutherland calculated the thermodynamics for SCM [13|] and found that the densities 
of occupied (pp) and unoccupied (pn) k satisfy the following relation: \pp{k) + puik) = 1. 
This is a statement of broken particle-hole symmetry for A ^ 1 and an analog of the Chern- 
Simons duality. pp{k) satisfies the following relation [|14]], (1 — App) A (l — (A — l)pp) A_1 = 



ppexp((e(k) — fi)/T) with e(k) = k 2 , which is identical to that of the ideal gas obeying Hal- 
dane's fractional exclusion statistics ||15|| . I emphasize here that the statistical distribution 



above is only for the pseudo-particles and is not related to the statistics of the real particles. 

I use three different names for the particles in the model — real, pseudo-, and quasi- 
particles. The real particles are, of course, the physical quantum particles described by the 
canonically conjugate coordinate and momentum variables {xj,pj}. The pseudo-particles 
are described by the pseudo-momentum operators (see Eq. fli~5|) ) whose eigenvalues are given 
by Eq. (|i~0|) . The quasi-particles are the elementary excitations of the system. Because the 
pseudo-particles form an ideal gas, the quasi-particles are essentially same as the pseudo- 
particles excited out of the condensate. The holes left behind in the pseudo-particle con- 
densate will be called quasi-holes. Hence, the name "pseudo" and "quasi" will be used 
interchangeably in some cases. 

If A = p/q, where p and q are relative primes, the motifs for eigenstates can be con- 
structed. The ones and zeroes in the motif mean occupied and unoccupied kj, respectively. 
From Eq. ( |TUD the following rules can be deduced for constructing the motif, (i) Total of 
N ones in the motif represent the pseudo-particles. Hence, the charge of pseudo- and real 
particle is same, (ii) The allowed number of zeroes between each pair of ones is p— 1 +nq, 
where n is an arbitrary non-negative integer. Of the p — 1 + nq allowed zeroes, p — 1 con- 
secutive zeroes are bound to each one while the rest of them are un-bound. The consecutive 
q un-bound zeroes represent a hole, (iii) A cutoff is introduced to make the number of un- 
bound zeroes finite. The natural unit for the pseudo-momenta is then 2w/qL which is an 
analog of the flux quantum. To give some examples, the ground state (a), allowed (6) and 
forbidden (c) states for A = 3/4 for N = 7 are represented, respectively, by 

(a) ■ • • 00000000100100100100100100100000000 • ■ ■ 

(b) ■■■ 00000000100100100100100100000010000 • • • 

(c) • ■ • 00000000100100100000100100100100000 ■ ■ ■ 

The motif (c) violates Rule (ii) and hence not allowed. The bound zeros which cause spacings 



between the ones can be viewed as flux attached to the pseudo-particles. 

The motif thus constructed is a powerful tool in visualizing the fractional statistics. 
Destroying a particle (or turning a one into a zero) creates p extra un-bound zeroes and 
in order to have integer number of holes in the condensate (Rule (ii)), minimum of q — 1 
extra quasi-particles must be excited out of the condensate leaving behind total of qp un- 
bound zeroes which break up into p holes. Hence, q particles leave behind p holes. This is 
a generalization of the Pauli exclusion principle which explains intuitively why the density- 
density correlation function has q quasi-particles and p quasi- holes when A = p/q. It will be 
shown explicitly that the only states that contribute to the DDDCF are indeed this minimal 
excitation. 

Before I begin to discuss the method used to calculate the DDDCF, some notations need 
to be defined. First, a diagram T>(k) is defined to be rows and columns of boxes labeled by 
'■ 1 < i < 1 < j < Ki}, where 1(k) denotes the number of non-zero Kj. The label 
i and j are row and column indices of the diagram with 1(k) rows of lengths Kj. Second, the 
conjugate of k denoted by k! is obtained from k by changing all the rows to columns in non- 
increasing order [IB]. Each row (column) corresponds to the quasi-particle (hole) excitations. 



Third, a generalization of factorial is defined by [a]£ = n(ij)g K ( a + (j ' ~ 1)M — (i — 1))- Using 
the above notations, the density-density correlation function at finite N and L is given by 

(o|p(mmo,o)[o> = A E M! _^m_ ^\ K \ x/ L e -,t E ^ 

where E K = Ef =1 k] + A Ef=i (N + l- 2j) Kj and j x K = ,,!/,•} - % + 1 + (/* - j)/X) (^ - 

i + («i — j + 1)/A). The product in [0']£ does not include the pair (z, j) = (1, 1). In deriving 
the finite size correlation function, I use the following two relations |]T7Hl8ll 

e^=;e?^(w). (12) 

i A | re [=n 

(k|k') a = j' A - -5 K (13) 

\ I /A ^[ jV+1 / A _ 1 ]A ^ ^ 

The first relation can be used to expand the density operator (1/L) Sj=i — x i) i n terms 
of the Jack symmetric polynomial. The second relation gives the normalization constants 
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of all the excited states parameterized by k. The DDDCF is now easy to calculate since the 
eigenstate k evolves in time only with a phase exp(— itE K ). 

For A = p/q, the coefficient [0']£ is zero if T>(k) consists of more than p columns or more 
than q rows. Hence, in the thermodynamic limit a local density operator acting on the 
ground state provokes only the minimal excitations consisting of q and p quasi-particles and 
holes. 

The large N expansion of Eq. flllD can be carried out at fixed density p = N/ L and the 
leading order term corresponding to the thermodynamic limit is given by Eq. (g). In this 
limit, a new super selection rule emerges and suppresses the states with \k{ — Kj\ ~ 0(1) 
to 0(1/N 2X ) or 0(l/N 2 ' x ) depending on the value of A. This means that the states with 
quasi-particles and quasi-holes with same momenta (velocities) are suppressed. There are 
some exotic exceptions to this rule. The details will be published elsewhere. 

The form of the ground state wave-function, Eq. fl8|), has lead Haldane to suggest that 
while the apparent statistics can be modified with a singular gauge transformation, the 
"natural" statistics of the SCM are fractional, and that the particle excitations carry charge 
1 and flux ir\ and the hole excitations charge — 1/A and flux —71 []nj. Indeed, if a singular 
gauge transformation [20[ is applied, the ground state wave function can be rewritten as 



*.=n£^to-*i*n*r ,( " M,/ '. <") 

j>l 1*3 Zl \ k 

where the apparent statistical parameter is now 9 = n/3. Unlike in the two-dimensional case, 
the transformed wave function remains as the ground state of the original Hamiltonian and 
so are all the other eigenstates. The DDDCF is also unchanged. The HPOGF, however, 
will be different for different choices of statistics (see, for example, 123). Hence, in order to 
calculate the HPOGF, it will be necessary to adapt Haldane's "natural" exchange statistics 
for the SCM particles. 

In order to consider the fractional exchange statistics, it is convenient to multiply the 
wave function \1> by an "ordering function" ip(xp 1 , . . . ,xp N ) which is just a book keeping 
device for the phase factor that depends on the permutation P and is set to unity for the 



fundamental region x± < X2 < ■ ■ ■ < xjv- (In the case of Fermions the function is just a 
product of Grassmann numbers.) Since the ordering function automatically keeps track of 
all the exchange phases, the particle exchange operator Pjj acting on the full wave function 
amounts to simply exchanging the indices i and j. If the "natural" statistics is chosen, the 
phases arising from ip and \I/ can always be set to cancel each other! 

In the original SCM there is no physical process which allows particles to exchange, (i.e. 
The wave function vanishes like I ) A few years ago, Polychronakos |22|| 

solved this problem by introducing an analog of the Yangian generator [|J 

Ti j = pj + in\/L'^2cot(7i(xi — Xj)/L)Pij, (15) 

where pj is the ordinary momentum operator, and showed that the Hamiltonian H = J2j ft] 
is fully integrable and is same as the SCM up to a modification A(A — Py) plus some trivial 
constant. The new operator ttj is the momentum operator for the pseudo-particles, and the 
momenta J2j kj corresponds to the eigenvalues of J2j ftj- This new Hamiltonian should be 
considered as the model of one-dimensional "anyons" with fractional exchange statistics. 

The single particle destruction operation on the ground state of iV + 1 "anyons" is, then, 
given by 

N 

9(x)\Q) N+1 = z- XN ' 2 H(z- Zj ) x zT x/2 \0) N , (16) 

where z = exp(i27rx / L) . Now, a similar technique used for the DDDCF can be employed 
to solve for the HPOGF. In this case, the contributing partitions have no more than p rows 
and q — 1 columns (i.e. p quasi-holes and q — 1 quasi-particles). Therefore, the "natural" 
exchange statistics of the real particles is fully compatible with the exclusion statistics of 
the elementary excitations. 

In the thermodynamic limit the one-particle Green's function (hole propagator) is given 

by 

(0|*tOM)*(0,0)|0) = pDe-^l[(rdx t ) fld'dyj) F(q - l,p, \\{x t , yj })e« Qx - Et \ 

i=i \ Jo ' j=i \ Jo ' 

(17) 

8 



where F(q — l,p, \\{xi,yj}) is given by Eq. (0) and D by 

A 2 f(?-i)rd + a) 

D - ^ {q - iy J' A(q - l ' p ' x) - (18) 

Q and E are same as before except for the number of Xj. At integer values of A (i.e. 
q — 1 case where only quasi- holes are excited), based on the equal-time results of Forrester 
23fl Haldane made a conjecture |1| which agrees with this formula. I conjecture that the 



minimal form factor for any two-point correlation function is given by F(m,n, \\{xi, t/j}) if 
the intermediate states involve only m quasi-particles and n quasi-holes. 

In conclusion, the SCM is shown to possess the fractional exclusion and exchange statis- 
tics. The motifs representing the full spectrum are constructed and used to demonstrate the 
exclusion statistics, explicitly. The fractional statistics in the SCM is also confirmed by cal- 
culating the exact dynamical density- density correlation function and one-particle Green's 
function (hole propagator) at any rational interaction coupling constant using the theory 
of Jack symmetric orthogonal polynomial. The details of the calculation including the full 
Green's function will be published elsewhere. 

While this Letter focuses on the fractional statistics aspect of the SCM, the method for 
calculating the correlation functions developed here could be of interest to a wide variety of 
people working on the disordered electronic system, the quantum chaos, the random matrix 
theory, 2D QCD, etc. 

I thank F.D.M. Haldane, F. Wilczek, Y. M. Cho, T. Hwa, C. Johnson, R. Narayanan, 
and R. Kamien for useful discussions. This work is supported by DOE grant #DE-FG02- 
90ER40542. 
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